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Abstract 

A previous paper constructed a kinematic basis for spin networks 
with planar or cylindrical symmetry and arbitrary polarization. This 
paper imposes a constraint which limits the gravitational wave to a 
single polarization. The spectrum of the constraint contains a phys- 
ically reasonable number of zero eigenvalues, and the zero eigenvec- 
tors can be constructed explicitly. Commutation of the constraint 
with the Hamiltonian is expected to lead to a further constraint. 
This new constraint is not investigated in this paper, but I argue it 
will be non-local, relating states at two or more neighboring vertices. 

PACS categories: 04.60.Pp, 04.30.-w 
I Introduction 



In a previous paper, I constructed a kinematic basis for a spin net- 
work with planar or cylindrical symmetry and studied the eigen- 
functions and eigenvalues of the volume operator [Tj. In that paper, 
the gravitational wave could have either polarization. In this paper 
I study the volume operator for a wave limited to a single polariza- 
tion. 

In classical gravity, and in quantum gravity based on local field 
theory, the assumption of a single polarization greatly simplifies the 
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discussion. If both polarizations are present the matrix of metric 
components (or triad components, if one is using connection triad 
variables) has a 1x1 subblock and a 2x2 subblock after gauge fix- 
ing. The 1x1 subblock contains the longitudinal component of the 
metric or triad (r component in the cylindrical case; z component 
in the planar case); the 2x2 subblock contains components trans- 
verse to the direction of the wave. In the single polarization case, 
the matrix simplifies further, becoming diagonal. In terms of ADM 
metric connection variables, the off-diagonal element g^y (or ge^)=Q. 
In the connection triad variables used in this paper, the equivalent 
statement is 

Tr(E y) = 0, (1) 

or the analogous statement for cylindrical coordinates. 

In local field theory, when the constraint eq. (Q) is commuted 
with the Hamiltonian, one obtains a further constraint which states 
that the conjugate ADM momentum vr^^ also vanishes. In triad 
connection language, 

= Tr(E^E'')Tr(EMJ) + (a;4->?/)-r^f; 

Y^y = -iTr((E ^ 9, E ^)E "). (2) 

The r term subtracts out the spin connection part of A. Imposing 
these constraints causes a degree of freedom and its conjugate mo- 
mentum to disappear entirely from the calculation, which simplifies 
considerably. 

In the spin network case, it is straightforward to impose eq. (jT)). 
This constraint resembles the volume operator, in that it is con- 
structed entirely from triads and is derivative free. If one makes 
the simplest assumption that only the field strengths (or other op- 
erators with derivatives) need to be non-local, then one can take 
the operator of eq. (H)) to be local. That is, the operator acts at a 
single vertex, merely reshuffling the holonomic basis states at that 
vertex. The constraint is no harder to handle than the volume op- 
erator, which is also local. (All E operators should be integrated 
over an area, but this is straightforward and I suppress the area 
integrations.) 

However, in a spin network approach, the constraint eq. is 
unacceptable. Connections A are not diffeomorphism invariant un- 
less integrated over an edge and included in a holonomy. It is also 
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much easier to construct a field theory based on holonomies than one 
based on A's, because holonomies may be described by a compact 
manifold, the Euler angles 0. 

Perhaps one could propose a spin network generalization of eq. 
that is local. (Promote the A operator to a holonomy which acts 
at the same vertex as the E operators in eq. 0) However, I will 
give two arguments that the holonomy (and perhaps some of the 
E as well) will have to be non-local. That is, the holonomy which 
regularizes the A in eq. will include segments joining two or more 
neighboring vertices. 

(First argument.) Since the constraint eq. (|21) is a momentum, 
the time derivative of the original constraint, the natural way to 
obtain (the spin network version of) eq. Q is to commute eq. ([Q) 
with the spin network Hamiltonian. It is not clear at present what 
to use for that Hamiltonian; but it must be non-local, in order for 
propagation to occur. Commuting even a local object with a non- 
local object is very likely to lead to a non local object. 

(Second argument). Consider the Killing vectors which define 
planar symmetry. (I am concentrating on the planar case. The 
discussion for the cylindrical case is virtually identical 0.) The 
space possesses two spacelike Killing vectors which commute, and 
one can choose coordinates so that these two vectors become d/dx 
and d/dy . This planar symmetry by itself does not imply one 
polarization only. Given only the existence of the two Killing vec- 
tors, all three transverse components of the metric, gxx, Qyy, dxy or 
all four transverse triads can be non-zero; this is enough to permit 
two polarizations. 

In classical general relativity, and in local quantum field theory 
treatments of gravity, one must demand hypersurface orthogonality 
of the Killing vectors in order to get one polarization . Hypersur- 
face orthogonality can be shown to imply that g^y and its canonical 
momentum vanish 

Hypersurface orthogonality is an intrinsically non-local idea. If 
a Killing vector field such as d/dx exists, then in each plane z = 
constant = zi there exist one- dimensional integral curves with the 
Killing vector as tangent vector. The plane at z = ^2 possesses 
similar integral curves, but no relation is implied between the two 
sets of curves at Zi and Z2. The additional assumption of hypersur- 
face orthogonality creates such a relation. The two sets of integral 
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curves must be normal to a common set of hyperplanes . The in- 
tegral curves can be described as curves of increasing x, and the 
hyperplanes as hyperplanes of constant x. 

Without hypersurface orthogonality, the integral curves of d/dx 
may twist in any manner as one goes from one value of z to the 
next. With hypersurface orthogonality, this twisting is strongly con- 
strained. 

In the spin network context, going from one value of z to the next 
is equivalent to going from one vertex to the next. Hypersurface 
orthogonality implies a relation between structures at two different 
vertices, necessarily a non-local relation. (End of second argument). 

If eq. (121) must be promoted to a non-local constraint, then within 
a spin network context, the one-polarization case is harder than the 
two-polarization case. In this paper I carry the discussion of the 
one-polarization case as far as I can, given that I have a kinematic 
basis but no Hamiltonian, therefore no way to construct eq. 0. I 
impose eq. ((T)) on the states of the kinematic basis, but not eq. (j2I). 
I then ask which states of the two-polarization case are ruled out by 
the constraint. 

For orientation, it is helpful to keep in mind an example from 
quantum mechanics, the two dimensional simple harmonic oscillator 
with Hamiltonian 

H = -(hV2m)[9V'9x2 + d^/dy^] + (k/2)(x2 + y^). (3) 

This problem is separable, with eigenf unctions | n^Uy > labelled by 
occupation numbers. Suppose I wish to impose the constraint that 
the oscillator possesses only a single "polarization", say it oscillates 
only along x. I can impose the constraint on the classical theory 
first, then quantize (constraint first). Alternatively, I can follow 
Dirac and quantize first, then impose the constraint (quantize first) 
For reasons to be explored in the next section, I will choose to 
follow the second, Dirac procedure and quantize first. 

Suppose I apply the Dirac procedure to the oscillator problem. I 
know the answer I want: the constraints should rule out all states 
except those of the form | n^iny = >. However, if I require 
physical states to obey 

y|>=Pyl>=o, (4) 

I find these requirements are too strong. From the theory of coherent 
states, there is no ket which satisfies even one of these equations 
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(Recall that both coordinate and momentum contain creation 
operators, which make it difficult for these operators to completely 
destroy any state.) 

I can try the milder constraint, 

<|y|>=o, (5) 

which is perhaps enough to guarantee a satisfactory classical limit. 
However, whereas the previous constraint was too strong, this con- 
straint is too weak. For example, it is satisfied by any state |> 
which has definite parity, including states with Uy non-zero. To get 
the answer I want, I must impose at least eq. (0) and 

<|y'l>=o. (6) 

Eq. gives the energy a satisfactory classical limit. I conclude a 
straightforward application of the constraint to a ket, the analog of 
eq. (j^, may not work; I may have to impose averaged constraints, 
eqs. ® and ©. 

For another example, this time from field theory, consider quan- 
tization of the electromagnetic field in Lorentz gauge. One must 
impose both of the constraints 

= <\d-A\>- 

= <| {d-Af |>, (7) 

in order to get the correct classical limit, including the correct clas- 
sical energy 0. Again, 

= (9-A|> (8) 

is too strong. 

These examples determine much of the discussion in section two. 
First I discuss the relative merits of constraint ffist vs. quantize 
ffist approaches, and I opt for the latter. I then impose the analog 
of eq. 0, 

Tr(E"Ey) |>=0. (9) 

In the oscillator and Lorentz gauge examples, constraints of this 
type have no solutions. However, in the spin network case I was 
able to find normalizable states satisfying eq. Q . The kernel of the 
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constraint is non-trivial. As a check, I impose the constraints in the 
average sense: 

= <| Tr(E y) \> 

= <\ [Tr(E "E y)]^ \> . (10) 

I find that the states which satisfy these constraints are the same as 
the states which satisfy eq. 0- 

II Dirac Quantization of the One-Polarization 

Case 

I can impose the constraint on the classical theory first, then quan- 
tize the single-polarization classical theory (constrain first); or I can 
quantize first, then impose the constraint as a condition on the quan- 
tized states (quantize first; Dirac quantization). If I constrain first, 
it is natural to follow Bojowald and introduce fields which have the 
constraint built in [7]. 

W^Sa = ^^[cos(a + /3)^i + sin(a + /5)^2]; 

E^^Sa = Ey[-sm{a + (3)Si + cos{a + (3)S2]; 

A^Sa = A,[cos{P)S, + sm{P)S2]; 

ApA = A,[-sin(/5)5i + cos(/3)^2]. (11) 

The first two lines, for example, relate the four "Cartesian" fields 
(E^ and E^) to three "polar" fields (the magnitudes E^,Ey, and the 
angle a + P). One degree of freedom is lost, but this is acceptable: 
the new polar fields have the constraint Tr(E ^'E ^) = built in. 
From eq. ()11|) the A's obey a similar constraint, 

TrA.Ay = 0, (12) 

which again removes one degree of freedom, and is easier to im- 
plement than the constraint of eq. 0- Bojowald carries out the 
canonical transformation from Cartesian to polar fields and finds 
the following pairs of canonical variables [7]. 

(E^cos(a) = py,Ay); 

{sma{A^P^ + AyPy) = P^,P); 

(E|,Af). (13) 
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turns out to be the Gauss Constraint. The angle a is not an 
independent degree of freedom, but is given by 

tana = P^/(AxP" + AyP^). 

These are classical field theory definitions. It is not clear how to or- 
der the operators in tan a after quantizing, or what the spin network 
definitions should be. 

The lack of precise definitions is already a problem at the kine- 
matic level, since the formula for the volume operator in polar co- 
ordinates involves a. 

= WzP'^py /{cosaf. (14) 

Since I was unable to supply any convincing spin network definitions 
of a or V^^, I determined to avoid polar coordinates as much as 
possible. Every time I attempted a "constrain first" quantization, 
however, I found myself introducing the polar variables, explicitly or 
implicitly. A constrain first quantization which avoids polar fields 
may be possible; but I was unable to find it. 

For the rest of this paper I will use the quantize first approach. 
This has the advantage that the first steps (gauge fixing, quantiza- 
tion, kinematic basis and dot product) have been worked out already 
in the previous paper; and the polar fields need not be used at any 
stage. 

I now impose the constraint, in its unaveraged form, eq. 0. I 
take the ket to be a linear combination of states in the kinematic 
basis. Note that the trace in eq. (jHl), is very similar in structure to 
the 2x2 transverse subblock of the (square of the) volume operator: 

= i(E^E^ -E^E^); 
^xy ._ 2Tr(E ^E y) 

= E^E^ + E^E^. (15) 

on the second lines I have used the linear combinations with simple 
U(l) transformation properties, 

E| = (E^ ± iEf.)/y2. 
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The two expansions in eq. (fT^ are identical, except for the extra i's 
and the extra minus sign in (V2)^. Therefore I can get the action of 
the constraint from a simple modification of the action of the (V2)^ 
operator. 

First recall the equation for the action of (V2)^. Write an eigen- 
function of this operator as a linear combination of states in the 
kinematic basis: 

I c(A; Lx, Ly, F) >= Sdc(D; A, L,, Ly, F)YL,^,YL,n.,. (16) 

A is the eigenvalue, c is the eigenvector with components {c(D;A, L^, Lj^,F)}, 
the Y's are the spherical harmonics constructed in pp, and 

F = {m^ + my)/2; 

D = {m,-my)/2. (17) 

The eigenfunction contains a sum over D only; the volume operator 
does not change F. From the previous paper, acting on this 

state leads to the following recurrence relation for the c's. 

2Ac(L>;A) = ig^ c{D - 1; \) - ig+ c{D + 1; \); 



g.{D) = J(L,-F-D + l)iL^ + F + D) 



\Ly + F-D + l){Ly-F + D)- 
g4D) = g-{D + l). (18) 

(For simplicity I have suppressed some of the arguments of the c's.) 
To obtain the action of E^'^, I use eq. ()15|) . I drop the i's and the 
minus sign in eq. (fTH|). Also, I replace the eigenvalue A by zero, since 

E'^y-^ = 0. 

= g_e(D-l)+g+e(D + l). (19) 

I denote the expansion coefficients by e, rather than c, because the 
c's (the volume eigenfunctions) in general will not satisfy eq. (jl9p . 

The following lemma is straightforward to prove. 

Lemma . If the eigenvector c(A) with components {c(D;A)} is a 
solution to eq. (fTHj). then the eigenvectors 

e(±A) = {(Ti)°c(D;A)} (20) 

are eigenvectors of eq. (fTUj) with eigenvalues ±A, i.e. 

± 2Ae(D; ±A) = g_e(D - 1; ±A) + g+e(D + 1; ±A). (21) 
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This lemma implies that (V2)^ and the constraint operator have 
exactly the same eigenvalue spectrum (although their eigenvectors 
are not the same). (For example, in the previous paper I proved 
that the non-zero eigenvalues of {¥2)"^ occur in pairs (A, —A); and 
the lemma indicates a similar pairing for the eigenvalues of the con- 
straint operator.) 

From the previous paper, {¥2)'^ has zero eigenvalues p. Hence 
the constraint operator has zero eigenvalues, and the constraint 
equation eq. Q has non-trivial solutions. 

From the discussion in the introduction, I should also investigate 
the averaged constraints, eq. (fTUj) . It is easy to see that there are 
a large number of solutions to the first, linear constraint: the com- 
binations I e(A) > ± I e(— A) > are sent into the orthogonal com- 
binations I e(A) > =F I e(— A) > by the constraint. (As at eq. (fTH|l . 
I e(A) > is shorthand for the eigenstate of the constraint operator 

However, the second, quadratic constraint in eq. ()10|) is more 
restrictive. Let the quadratic constraint act on an arbitrary state 

|>=^a(A)|e(A)>. 

A 

The expectation value of the quadratic constraint is then 

SaA^ I a(A) p . 

This is not zero unless all of the A are zero. The averaged constraints, 
eq. (fnH). give the same result as the unaveraged constraint, eq. 

In the oscillator example discussed in the introduction, the im- 
position of a constraint lowers the dimension of the Hilbert space 
signifigantly. The unconstrained oscillator states are labeled by two 
integers {nx,ny), and there are (countable infinity)^ states in the 
Hilbert space. After the constraint restricts the states to = 
only, the number of admissible states drops to (countable infinity)^. 
Roughly speaking, this is the square root of the original number of 
states. 

Something similar happens in the present case. Consider the set 
of kinematic basis states {YL^rn^YL^my^Wrrix, my}. Before imposition 
of the constraint, there are 

(2Lx + l)(2Ly + 1) unconstrained states 



9 



in this set. After imposition of the constraint, the set contains only 
those hnear combinations which are eigenstates of the constraint 
with eigenvector zero. From the lemma, the volume and constraint 
operators have exactly the same number of zero eigenvalues. The 
necessary condition for the volume operator to have a zero was 
worked out in the previous paper [P. The eigenvectors are sums 
of the form eq. sums over D with F held fixed. There will be a 
zero eigenvector iff the number of D's in the sum is odd. Therefore 
I can compute the number of zero eigenvectors by constructing the 
following rectangular lattice of points. Make the x-axis rrix and y- 
axis rriy] draw one point with coordinates (m^, my), for each allowed 
pair of m's, (2La. + l){2Ly + 1) points in all. 

In this lattice, diagonals at 135 degrees are lines of constant F = 
{rrix + my)/2, while diagonals at 45 degrees are lines of constant D 
= {rrix — my)/ 2. The number of zeros is the number of 135 degree 
diagonals which contain an odd number of points (corresponding 
to an odd number of D's in the sum over D with F held fixed). 
After some experimenting with specific examples, one arrives at the 
general formula for the number of diagonals with an odd number of 
points = number of states surviving the constraint. 



We have gone from order (2L)^ states to order (2L)^ states. This is 
roughly a square root, similar to the oscillator example and therefore 
not unreasonable. 

In the previous paper I computed the zero eigenvectors of the 
volume operator pp. Given that result, and the lemma, it is pos- 
sible to compute the zero eigenvectors of the constraint. From the 
previous paper, the zero eigenvectors of the volume operator have 
components 



for D = max D, max D - 2, max D - 4, . . ., min D; and c(D) 
= otherwise. From the lemma, eq. ()2()|1 . the components of the 



2 max(Lx, Ly) + 1 constrained states 
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corresponding zero eigenvector of the constraint are just 

e(D; A = 0) = (i)°c(D; A = 0). 

(The two states e(±; A = 0) in eq. give the same zero eigenvec- 
tor, except perhaps for an overall minus sign, because c(D;A= 0) is 
zero for every other value of D.) 

These zero eigenvectors illustrate a point made earlier, that vol- 
ume and constraint operators have (the same eigenvalues but) very 
different eigenvectors. Zero constraint eigenvalue does not imply 
zero volume eigenvalue. An eigenvector of the constraint opera- 
tor with zero eigenvalue is a linear combination of eigenvectors of 
the volume operator, with eigenvalues which are in general different 
from zero. Although the volume operator and constraint have the 
same eigenvalue spectrum, they do not commute. 

I return to the point made in the introduction: in classical or 
quantum field theory, the single polarization case is simpler; whereas 
for spin networks, the single polarization case is harder. For exam- 
ple, experience with field theory might suggest using eigenf unctions 
of the constraint as a basis; however, these eigenf unctions are not 
eigenfunctions of the volume, and the spin network Hamiltonian 
would be hard to evaluate in that basis. 

There is a further complication. The constraint I have been 
studying in this paper is a local constraint. The Hamiltonian acting 
on an eigenfunction of this local constraint will produce another local 
constraint eigenfunction only if the state satisfies the spin network 
version of eq. Q. As argued in the introduction, that constraint 
is likely to be non-local. Given the problems with choosing a basis 
and the existence of a non-local constraint, it would not be surpris- 
ing if solutions for a general polarization were obtained first, before 
solutions for a single polarization. 

There may be ways around these complications, however. In 
the classical limit, the spin vectors associated with the and Ay 
holonomies will have sharp expectation values with minimal uncer- 
tainties. The local constraint implies that the two spin vectors are 
perpendicular; while the volume eigenvalue is given by the cross 
product of the two spin vectors. It is not possible for the classical 
state to be simultaneously an exact eigenfunction of both the vol- 
ume operator and the local constraint; but in the classical limit it 
should be possible to find states which are eigenfunctions of both 
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constraints to very good approximation. Further, the non-local con- 
straint, when combined with the Hamiltonian, may simplify the 
Hamiltonian. It is too early to tell. Much work remains to be 
done. 
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